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Abstract 

We consider alloy type random Schrodinger operators on a cubic lattice 
whose randomness is generated by the sign-indefinite single-site potential. 
We derive Anderson localization for this class of models in the Lifshitz tails 
regime, i.e. when the coupling parameter A is small, for the energies E < 
-CX 2 . 



1 Introduction and main results 
1.1 Introduction 

The prototypical model for the study of localization properties of quantum states 
of single electrons in disordered solids is the Anderson Hamiltonian Ha on the lat- 
tice Z d . It consists of the sum of the finite difference Laplacian that describes the 
perfect crystal and a multiplication operator by a sequence of independent identi- 
cally distributed random variables that emulates the effect of disorder. The basic 
phenomenon, named Anderson localization after the physicist P. W. Anderson, 
is that disorder can cause localization of electron states, which manifests itself 
in time evolution (non- spreading of wave packets), (vanishing of) conductivity 
in response to electric field, Hall currents in the presence of both magnetic and 
electric field, and statistics of the spacing between nearby energy levels. The first 
property implies spectral localization, i.e. the spectral measure of is almost 
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surely pure point in some energy interval, and almost sure exponential decay of 
eigenfunctions there. 

In this paper we study the properties of the more general class of systems, 
called the random alloy type model, in the three-dimensional setting. The action 
of the corresponding Hamiltonian H^ on the vector y/ e / 2 (Z 3 ) is described in 
equation below: 

(tfJvOW : = -^(A ¥ )(x)+XV C0 (x)xjf(x). (1) 
Here A denotes the discrete Laplace operator, 

(AvO(x) = £ V(x+e) - 6yf(x), (2) 

e<El?, |e[=l 

and V a stands for a random multiplication operator of the form 

Vco(x) = £ COiu(x — i). (3) 

iekl? 

Here fcZ 3 denotes the set of all points i on 1? which are of the form 

i = {hziMziMz-i) 

for an arbitrary vector z=(zi, zi, Z3) 6Z 3 and a given vector k=(k\, k%, £3) G I? . 
The parameter A conveniently describes the strength of disorder. We note that 
the standard Anderson Hamiltonian Ha corresponds to the choices u(x) = 8 X ; k = 
(1,1,1) (so that kZ 3 = I?). Here 8 X stands for for the Kronecker delta function. 

The localization properties are known to hold for Ha in each of the following 
cases: 1) high disorder (A 3> 1), 2) extreme energies, 3) weak disorder away 
from the spectrum of the unperturbed operator, and 4) one dimension. Most of 
the mathematical results on localization for operators with random potential in 
dimensions d > 1 have been derived using the multi-scale analysis (MSA) initiated 
by Frohlich and Spencer ifTOll and by the fractional moment method (FMM) of 
Aizenman and Molchanov [0. The Anderson localization problem for Ha has 
been studied comprehensively, see for example [fT6l and references therein. 

For a general single site function u the situation is more complicated. If all 
the coefficients u{x) have the same sign, the dependence of the spectrum on V is 
monotone - property that obviously holds true for Ha- Localization in such sys- 
tems (in all regimes above) is relatively well understood by now, using the meth- 
ods developed for the original Anderson model. There is however no physically 



2 



compelling reason for a random tight binding alloy model to be monotone, and the 
natural question is whether Anderson localization still holds in the aforementioned 
regimes in the non-monotone case, i.e. when u(x) are not all of the same sign. 
Mathematically, the problem becomes especially acute when u := Y# u ( x ) = 0> 
see the discussion at the end of this section. In the strong disorder regime the 
recent preprint extended FMM technique for a class of matrix Hamiltonians 
that include//^. 

In this paper we study the localization properties of the alloy type models in 
the so called Lifshitz tails regime, namely localization of the states with energies 
that lie below the spectrum of the free Laplacian, for A. <C 1 . The occurrence of 
localization for Ha at energies near the band edges at weak disorder is related 
to the rarefaction of low eigenvalues, and was already discussed in the physical 
literature by I. M. Lifshitz in 1964, see Section 3 in [12111 . and ll22l . 

As far as the rigorous results are concerned, there is an extensive literature 
devoted to the proof of Lifshitz tails as well as Anderson localization in this regime 
for the original Anderson model Ha- In this case 'mfa(HA) = CX almost surely, 
where the constant C is smaller than zero provided that the random variables (O x 
are i.i.d. and assume negative values with non zero probability (we will specify 
the assumptions on the randomness later on). One then is interested in showing 
that there is a non trivial interval / of localization at the bottom of the spectrum 
of Ha, with / = [CX,Eq(X)]. Among results in this direction, let us mention the 
work of M. Aizenman (TJ that established localization for Eq(X) = CX +0(X 5 ^ 4 ). 
This was later improved by W-M. Wang ||27ll to Eq(X) = CX + O(X) and then by 
F. Klopp OH to Eq(X)= CX 1 / 6 , with C < 0. Motivated by the unpublished note 
of T. Spencer 11261 . one of the authors flU derived the localization for the interval 
/ characterized by Eq(X) = CX 2 , with C < (it is expected on physical grounds 
that beyond the spectrum of Ha consists of delocalized states slightly above the 
threshold established in the latter paper). The proof utilizes the diagrammatic 
technique also employed in the current work. In fact a large portion of the effort 
spend here involves reducing the problem to the one studied earlier in [6], and 
generalizing the techniques presented there. 

Less is known in the case of a general single site potential. When the average 
u of the single site potential is not equal to zero, mfo(H^) < CX almost surely, 
with C < (see Section 5.1 of lfT8l ). F. Klopp established the region of localization 
with Eq(X) = CX 1 / 6 in this case, [fT9l . We show here that it persists at the energy 
range E < CX 2 regardless of the value of u. However, for u = the spectrum 
starts at CX 2 , and in order for the result to be non trivial in this situation, we have 
to show that C < C. At the end of this section we will describe the dipole model 
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for which such condition can be proven. 



1.2 Assumptions 

1.2.1 Randomness 

(gf) The values of the random potential {co ; } are i.i.d. variables, with even, 
compactly supported on an interval J, and bounded probability density p. 
We will further assume (without loss of generality) that J is centered around 
the origin and that the second moment satisfies E cof = 1 . The function p is 
a-Holder continuous: 

\p(x)-p(y)\ < K\x-y\ a max(lj(x),lj(y)), (4) 
with a > and where 1/ stands for a characteristic function of the set /. 

1.2.2 Single site potential 

We will focus our attention on two somewhat extreme cases: 

( G) Overlapping setup: The vector k is of the form k = (1 , 1 , 1 ), so that kl? = 
1?. This case corresponds in some sense to the maximal random setting, 
and here we will impose rather mild conditions on the single site potential 
function u. Namely, we will assume that in this case u(x) decays exponen- 
tially fast: 

< Ce- A ^. (5) 

{JV) Non overlapping setup: The vector k is such that {© — i} n © = for all 
kZ 3 , and supp u =: is compact. This setting correspond to non 
overlapping random potential. We will denote the corresponding primitive 
cell ©, i.e. © C © and the translates of by i e kl/ tile 7Lr . 

1.3 Notation and quantities of interest 

Let e(p) denote the dispersion law, associated with the Fourier transform of the 
Laplacian, (^Af)(p) = —2e(p)f(p), where 

f(p):=W)(p):= I^e-^fin), p e T 3 := [-1/2, 1/2] 3 , 
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with its inverse 

g(n)= f d 3 pe i2 ^ n f(p). 
Jt 3 

One then computes 

3 

e(p) = 2 £ sin 2 (^p-e a ), (6) 

where e a is a unit vector in the a direction. The spectrum of the unperturbed 
operator is absolutely continuous and consists of the interval [0, 6] . 

In what follows we will denote by A(x,y) the kernel of the linear operator A 
acting on / 2 (Z 3 ) (that is A(x,y) = (8y,A8 x ) = (y\A\x), where 8 X is an indicator 
function of the site x El?, and (•, ■) denotes the inner product of / 2 (Z 3 )). We will 
use the concise notation / in place of ff-py whenever it is clear from the context. 

The paper is devoted to the investigation of the properties of for a typical 
configuration CO in a weak disorder regime, namely at the energy range E < Eq, 
where 

1. For Case 0, 

E = Eg = -2A 2 ||«||i-2A 4 ||i2||i, (7) 
for X > being sufficiently small. 

2. For Case jV , 

E Q = Ef = -4nA 2 || M ||i{(6-2£) diam ®|0| + l}, (8) 

where is a primitive cell described in Assumption (jV), and |0| stands 
for the cardinality of the set 0. 



Diagramatic expansion and self energy. The quantity of the most interest is 
the typical asymptotic behavior of the so called Green function (also known as the 
two point correlation function, the propagator) 

R E+ ie(x,y) = (H^-E-ie)- l (x,y) (9) 

in the limit e \ 0. It plays a crucial role in determining, for instance, the conduc- 
tivity properties of the physical sample (whether it is an insulator or a conductor 
at a given energy band). On a mathematical level, investigation of the propagator 
can yield an insight on the typical spectrum of at the vicinity of E. 
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The technical assertions in this paper are proven using Feynman diagrammatic 
expansion for Re+ie around the unperturbed resolvent, i.e. the one with X = 
(Section [5]). The rate of convergence for this expansion in the limit e — > 
depends strongly on the value of E, and sets our limitations for the length of the 
interval / where we can prove localization. One can eliminate certain terms in this 
expansion (the so called tapole contributions) that are especially problematic from 
the convergence point of view, by modifying the unperturbed Hamiltonian. The 
corresponding addend a is called the self energy of the model. 

In the case (0) we define the self energy by the solution of the self-consistent 
equation 

a(p,E + ie) = X 2 f d\ , r J g( f ~ g) f -. (10) 

Ky ' .h e(q)-E-i£-a(q,E + i£) 

The relevant properties of the solution of (fTOb are collected in Appendix [Bj In 
particular, it has a single valued solution a(p,E + ie) for all e, all p G T 3 and all 
values of E that meet the condition ©. The function o satisfies 

Hallo. < min(-E -2X 4 \\u\\i,2X 2 \\u\\l). 

In the case {^V), the introduction of the self energy term requires some addi- 
tional preparation. Let n = \fa\. Let {*;}" = i be some enumeration of the sites of 
the primitive cell 0. Let D be the diagonal n x n matrix with Da = u(xi), and let 
@ be its periodic extension to £2(1?), i.e. 

@(x,y) = u(x mod fa) 8 x - y . (11) 

For any n x n matrix a, we construct the periodic operator E acting on Z 3 by 
defining its kernel as 



Oij x,y G + / for some I e £Z , x mod0 = ^, y mod& = Xj 
otherwise 

(12) 

We can now define annxn matrix S given by 



(-A/2- E-ie-L) 1 xj); Xi,Xje®. (13) 



Then the self energy term o in the case {JV) is going to be the nxn matrix, which 
satisfies 

o = X 2 DSD. (14) 
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The solution of (fl4l) enjoys properties similar to the one of (flOl ), namely it is 
unique for all E < Eq (where Eq is given by ©) and |e| < 2n X 2 ||w||L and satisfies 

|| cr || < 2nX 2 \\u\\i. (15) 

We defer further discussion of the properties of o to Appendix|B] 

1.4 Main result 

The hallmark of localization is a rapid decay of the Green function at energies in 
the spectrum of H a , for the typical configuration co. This behavior can be linked to 
the non-spreading of wave packets supported in the corresponding energy regimes 
and various other manifestations of localization. Our main result, Theorem \T\ 
below, establishes this behavior of the Green function at the band edges of the 
spectrum, by comparing it with the asymptotics of the free Green function. 

Theorem 1 (Anderson Localization for Lifshitz tails regime). For as above 
that satisfies Assumption (srf) and either {(J) or for any v > there exists 
Ao(v) such that for all X < Xq(v) the spectrum of H a within the set E <Eq — X 4 ~ v 
is almost-surely of the pure-point type, and the corresponding eigenfunctions are 
exponentially localized. 

1.5 Discussion 

It should be noted that our method works most effectively when the average u 
of the single site potential is not equal to zero. In this case infa(//^) < CX 
almost surely, with C < (see Section 5.1 of [18]). Hence for X sufficiently 
small Theorem Q] establishes the localization at the bottom of the spectrum of 

. We are only aware of one result on the Anderson localization in the regime 
discussed here for a non sign definite single site potential: F. Klopp proved the 
weak disorder localization for E < —A 7 / 6 in three dimensions, |fT9ll . Since for 
u = case — inf o(H^) = 0(X 2 ), his result does not provide an answer on whether 
the bottom of the spectrum is localized or not. 

When the average of the single site potential vanishes, we expect our method 
to yield the non trivial result when the minimizing configurations of the random 
potential look flat. That is essentially a reason why we can cover say the dipole 
potential in the {&) case below, the observation we owe to Giinter Stolz. In this 
case the expansion around the free Green function is a sensible procedure to do. 
However, in the case (^V) the sufficient symmetry of the single site potential can 
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cause the minimizing configuration to become periodic. For reflection symmetric 
u this was shown in [41 (for the continuum analogue of H^). In this situation, the 
free Green function does not capture the right features of the problem, and as a 
result our method fails to achieve the non trivial result in this case. It is worth 
noticing that exploiting the specific knowledge about the minimizing potential, 
one can show Lifshitz tails and consequently Anderson localization for such re- 
flection symmetric single site potential [20]. 

In this paper we consider the cubic lattice, that is d = 3 case. Similar (in fact 
better) results can be established for a higher dimensional case, but not for d < 3. 
Mathematically it is related to the nature of the point singularity of the propagator 
e(p) at zero energy - it is integrable for d > 3. This fact allows us to control the 
underlying Feynman series. 

We end the discussion with an example pertaining to the case w(0) = where 
the technique developed in this paper allows to get a meaningful result by improv- 
ing the bound on the threshold energy Eq\ 

Consider the single site potential u d of the dipole type, i.e. 

f 1 x = 

Ud(x) = < -1 x = e\ . (16) 
[0 otherwise 

Proposition 1. For the single site potential u d we have 

inf<7(f/£) < -2X 2 + 0(X 3 ) 

almost surely. The statement ofTheorem\l\holds true for all energies E that satisfy 

E<E d := -(1+A)A 2 . 

2 Outline of the proof 

We will establish Theorem \T\ using the multiscale analysis (MSA) method. It 
requires two inputs: The initial volume and Wegner estimates. 

2.1 Diagrammatic expansion 

The first ingredient in the proof of Theorem [Dean be established in the framework 
of Feynman graphs perturbation technique. Throughout the text various Green 
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functions will appear. We will denote by G E (x,y) the free Green function, i.e. 

G E (x,y) := (x\ (-A/2-E)- 1 \y) . (17) 

We characterize its relevant properties in Appendix [A] It will be used in some of 
the proofs as comparison with the full Green function R E +i e (x,y), defined in ©. 
Whenever it is clear from the context, we will suppress the energy dependence of 
Re+ie, and just use R (respectively R(x,y)) for the full resolvent (the full Green 
function). 

The following representations for a Green function R(x,y) provide the key 
technical tool for us: 

Lemma 1 (Decomposition of R(x,y) in (jV) case). Let E^ = —E + Eq/2, and 
let 



8,^ := y/E*r/3. (18) 

Then for any integer N and energies E < Eq with Eq satisfying ([8]) we have the 
decomposition 

N-l 

R(x,y) = £A„(*,y) + £ A N (x,z)R{z,y) , (19) 

n=0 zG Z 3 

with Ao(x,y) = R r (x,y) (the latter quantity is defined in (|52l) below), and where 
the kernels A n , satisfy bounds 

E\A n (x,y)\ 2 < (4n)!Z?> fc(£» e' 8 ^^ , n> I; (20) 

/ J 2 \ N / 2 

E\A N (x,y)\ < y^ANV. f C(£» -=\ e'^^ 2 , N> 1; (21) 

where C(E* / ^) = K |(a>| \\@\\ in 9 (E^) and K is some generic constant. 
The zero order contribution Aq satisfies 



\A (x,y)\ < G E *Jx,y) < Ce [x max , (1 + l^-yl)- 1 ) (22) 

for all x,y G Z 3 . 
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We now formulate the parallel result for the case {&). To this end, we intro- 
duce some additional notation first. For the parameter E* ff that satisfies — E +Eq > 
E^>0 with £,f defined in ©, we set 



~pff p 174 



(23) 



V6n 

Lemma 2 (Decomposition of R(x,y) in {iff) case). For any integer N and energies 
E < Eq with Eq satisfying © we have the decomposition 

N-l 

R{x,y) = £A n (x,j) + £ A N (x,z)R(z t y) , (24) 

n=0 ze z 3 

with Ao(x,y) = R r (x,y) (the latter kernel is defined in (|65l ) below), and where the 
(real valued) kernels A n , satisfy bounds 

E|A„(^)| 2 < (4n)\E* ff (c(E* ff ) -^=\ e-W x -y\, n>l; (25) 

E \A N (x,y) I < v^V)! \ C(E* ff ) e -**\*-y[/2 jN>1; (26 ) 

where 

C(E* ff ) = K(\\u\\^+A- 4 d ff )\n 9 (E* ff ) 

for some generic constant K and A being a parameter introduced in ©. 
The zero order contribution Aq satisfies 

\A (x,y)\ < 2e^ lx - yl (27) 

for all x,y G Z 3 . 

Remark 2.1. There is a certain balance between the parameters E* ff and 8$ in 
the above assertion. Namely, increasing the former improves on the pref actor in 
front of the exponential decay in (|25l) . but since it decreases the latter, the control 
over the rate of the exponential decay becomes poorer. 

Remark 2.2. The representations (fT9b and (|24|) are resolvent type expansions ( see 
Lemma^below for details). If one applies the rough norm bound on each factor 
of the resolvent there, the denominator in (|20l) and (|25l) will contain E* rather 
than its square root. The improvement is achieved using the Feynman diagramatic 
technique ( Section [3]). 
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One then looks for the optimal value ./V to stop the corresponding expansion - 
note that the increasing factor of (4iV) ! in Afj(x,y) competes with the decreasing 
factor (X 4 E*) N / 2 . 

The choice E* = A 4 ~ v /2 has the effect that 

A 2 

C{E*) -j= < X Bv , 0<B<1, (28) 

which suffices to control - (HB (respectively ([25]) - d2"gTT). We note that 
in the range of energies E < Eq — A 4 ~ v , the above choice for E* ff implies 8 > 
A 4 ~ v / {2k) 2 . It turns out that the appropriate choice for N should satisfy 

(see the next section for details). In terms of the A - dependence, it corresponds 
toN~ X~ bv for b<B. 



2.2 Wegner estimate 

The initial volume estimate that enters into MSA requires us to get rid of the 
imaginary part £ of the energy, present in the formulation of Lemma [2l To this 
end, we will use the Wegner estimate below, that itself is an important ingredient 
of MSA. It will be established using the idea of F. Klopp ifTTl . While it has the 
correct (linear) dependence on the length of the interval /, the dependence on the 
volume of the estimate below is not optimal. The optimal, linear dependence on 
A for the continuum models with absolute continuous density p was developed 
in |fT2|. but the trade-in is that the /-dependence in the Wegner estimate thereof 
is worse than ours. The derivation below has an advantage of being completely 
elementary and the estimate itself is sufficient for our purposes. 

Theorem 2 (Wegner estimate for Holder continuous densities). Let I be an open 
interval of energies such that 

Di := dist(/,a(-A/2)) > 0. 

Then for either O or ,jV cases we have 

ETrP/(#£' A ) < C\l\\t^{DiY x , (29) 

where H m denotes a natural restriction of to A C Z , the constant oc is 
defined in (&/), and the constant C depends on J and K. 
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Remark 2.3. The dependence on volume A blows up as Of \ 0. Thus the tech- 
nique doesn 't work when p is concentrated in finitely many points. See Bourgain 
and Kenig ^ for the Bernoulli alloy type model, where p is concentrated on 
{0, 1}, in continuous settings. 

The rest of the paper is organized as follows: We derive the main results based 
on the technical statements above in Section [3l We then describe the procedure 
which allows us to have a tighter control over the resolvent expansion in Section 
|4l We prove Lemmas Q] and [2] in Section \5\ Auxillarly technical statements are 
collected in Appendices lAl and iBl 

3 Proofs of the main results 

Proof of Theorem [7] 

The proofs of the auxiliary statements, namely Theorem [2] (accordingly Lemmas 
CD©, will be postponed until later in this section (until Section[5]). 

We prove the assertion simultaneously for {JV ) and {&) cases, so we are going to 
drop subscripts .jV or G until the very end of the proof. Let us denote by the 
natural restriction of H x to A C Z 3 , namely, (i, j) = H x (z, j) if (z, j) G A x A 

and Ha l (i,j) = otherwise. Let A c := Z 3 \ A, and let dA = {i\ 3 j's.t. (ij) G 
A x A c , dist(ij) = 1} be the boundary of the set A. We define the decoupled 
Hamiltonian H\ to be 

H K = H% X ®H^\ 

and will denote by R a (E) the corresponding resolvent (i.e. Ra(E) = (H^ — E) 
For L > and x G Z we denote by x = {y G Z d : \x — y |oo < L} the cube of side 
length 2L. Our first objective is to derive the bound for E\Ra l JE + ie;x,w)\, for 
w G dA. To this end, we observe 

E\R ALiX (E + ie;x,w)\ 

< E\R(E + is;x,w)\ + E\R Alx (E + i£;x,w) - R(E + ie;x,w)\ 
= E\R(E + ie;x,w) \ + E\ (R(H a - H Alx )R Alx ) (x, w) \ 
< E\R(E + ie;x,w)\ 
+ E £ {RiE + ieixtflliRj^iHa-H^fcwy. (30) 
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We can estimate 



I \(RA L ^(Hco-H AL ^)(k,w)\ < %dA\ = C=-, (31) 



hence 



E|tf A (£ + ie;x,w)| < C— max E\R(E + ie;x,k)\ . (32) 

£ kedA c Lx 

On the other hand, for any & 6 dA c Lx , we have dist(x,fc) = L+ 1, and Lemma Q] 
(respectively Lemma [2]) ensures that 



2V-1 

E|tf(E + /e;jc,*)| < £E|A„(jc,A:)| + Y, E \^N(x,z)R(E + ie;z,k)\ 



n=0 
N-l 



N-l / 
n=0 \ 



A 



2 \«/2 



,-«k-*l/2 



+ I 

zez 3 



A 



\/E* 

2 \^/2 



,-«|*- Z |/2 



JV— 1 / 

< e - SL/2 £ y/E*(4n)l C(E* 



A 



2 \"/2 



+ c 



2 \*f/2 



(33) 



Choosing 



(ANY 



C(£*)A 2 



(34) 



one obtains, using the Stirling's approximation, that the summation over the index 
n is bounded by a constant and 



(4N)\ 



C(£*)A 2 N N 
\fE* 
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-N 



Hence, for such a value of N we have 



E\R(E + ie;x,k)\ < C ( e ' 6L/2 + J^j 



Combining this bound with (|32l) . we obtain 



E\R A ^(E + i£;x,w)\ < C — 



e SL/2 + e 



-Nl 



£8 



U*I=[E-e l /*,E + e l /*],and let 

G(7) := {coe& : o{H Al J HI = 0} . 
For any co E G{I) we have by the first resolvent identity 

| R Alx (E + ie;x,w)- R Al ^ {E;x,w)\<£ 1 ' 2 . 
Pairing this bound with (1361) and using Chebyshev's inequality, we get that 



Prob I (D e G(7) : \R Alx (E;x,w) \ > C-^ 



r 8L/2 + e 



-N 



£8 



ProbUR ALx (E;x,w)\>C-^ 



e SL/2 e 



-Nl 



£8 



+ £ 



1/4 



(35) 



(36) 



(37) 



+ £ 1 / 4 

The Wegner estimate (l29l implies that 

Prob{cr(// AL Jn/^0} < C\I\\A L , X \^(D,)- 1 = Ce 1 / 4 ^/)- 1 /. 3 ^ . (38) 
Combining (1371 ) and (1381 ) we arrive at 



< Ce 1 / 4 ^;)-^ 3 ^. (39) 



We are now in position to set the values for the various parameters in the above 
formula, in terms of the single parameter A . We first note that since E < Eq — 
A 4 ~ v with < v < 1 and E defined in © - ®, it is allowed to choose E* = 
X 4 - v /2. It then follows from ([28]) and (El) that N ~ A~ B ' V with < B' < 1/4, for 
X small enough. Next, the parameter 8 originating from Lemmas[T|and[2]satisfies 
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5 ~ A 2 ~ v / 2 for A small enough. We finally choose L = A 2 and e = e~ 5 ^ Bv 2 _ 
Plugging it all into (l39l) we obtain that for the small values of A the following 
initial volume estimate holds true: 

C i—B'v/2} i — fl'v/2 

Prob | |fl Aji _ 2 {E;x,w)\ >e~ x \ < e~ x . (40) 

The initial volume estimate (|40l) together with the Wegner estimate (|38l) provide 
the necessary input for MSA for small A, and the result follows from say Theorem 

2.4 of mo. 

□ 

Proof of Theorem |2] 

In the sequel we will use the fact that any a-Holder continuous non negative 
function p admits a Lipschitz approximation by means of a non negative function 
pK such that px is AT-Lipschitz supported on J, and 

IIP— Px-lloo < CKT^. (41) 
It is then follows from ||p || i = 1 that 

||Pa:||i < 1+C(K^ \J\. (42) 

Observe that for any random quantity F a that depends on the random variables 
COi, with i E A, we have for any 8 

EF ro = f FcoHpi^dcOi = ^ [ 1+d v^dv [ F^YlP^dcbi. (43) 

J i€A 011 J ieA 

So, in order to evaluate E Pj(H^), we can first integrate over the fictitious ran- 
dom variable v. We will rely on the following simple statement. 

Lemma 3. Let A, 5 be hermitian nx n matrices. Let {A;t(A)}^ =1 , respectively 
{A^(5)}^ =1 be the set of corresponding eigenvalues in the ascending order, and 
suppose that 

0<a = Ai(£)<A„(£)=/3. 

Let I,J be the intervals [a,b] and [c,d] C R+ accordingly and let Pj denote the 
characteristic function of the interval I. Then we have 

jTrP!(A + xB)dx < a 1 \I\ TrP ; (A) , (44) 
where I = [a — [5d, b — ac\. 



15 



For all v G [1,1 + 5] we can use (I4TT) to bound 



piy&i) < p K (y(bi) +CKT^lj(va>j) 

< pK(cbi) + (cK^e + K8 \J\ ) If (obi) 



(45) 



where J' = (1 + 8) J. 

If E is a middle point of /, we can write Pj(H^) = Pj (H^ —E), where Iq is 
centered at origin and has the same width as /. Let B := —^A — E, then —E = 
B + vXv£. We note that B satisfies the same properties as its counterpart in Lemma 
[3] above, with a = Dj (from Lemma |2]) and /3 = 6 + Dj. Since 



P k (B + vXv£) = P v - % (v- l B + Xv£) < P h) (v- l B + XV l 
for v G [1, 1 + 8], we can estimate 

•1+5 

f(6Di)d(bi I 



EP T (H 



0) , 



< i [HfitydGH [ 1+ v^P^B + XV^dv (46) 
J ieA Jl 

< I (1 + S) |A| Jl\f(dH)ddH fj P h (v- l B+XV%)dv 
= ^(1 + 5)I A I Jjjfityd&i £ +3 P k (xB + Xv£)x- 2 dx 

< I(l + 5) 2+ l A l fHf{dk)ddH f Pj^xB + XV^dx, 



where 



J s = [(1 + 5)- 1 ,! 
Applying Lemma [3] with the choice A = X V A , we get 

EP/(4 A ) < U\ + 8) 2+ ^(D!)- l \I\ /n/(^)^TrP /5 (Ay A ), (47) 



where 



|/| 

Y 



-(6+D/),W-D7(l + 5)- 



(48) 
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Since TrP /g (AV A ) < |A|, we can use (SB and (g2j) to bound 

EP/(/4' A ) < i (1 + 5)2+1^ (Z> z )-i |/| |A| [nm)d<&i 

J i'eA 

< ^(l + 5) 2+2 l A l (l + (2C^+^|7|)|7|)' A| (Dj) 



-l 





|A| 


\J\ 


|A| 



< 



8D] 

where in the second step we have used (l45l) . Choosing 
K = (\J\\A\)^ ; 8- 1 

we get the desired bound (|29l ). 
Proof of Lemma \3\ 



1 + A \~\J \~ A 



□ 



The result follows from two consecutive applications of Weyl's theorem, which 
states that for any pair of n x n hermitian matrices C and D we have 

X k {C)+Xx{D) < X k (C + D) < X k (C)+X n (D). (49) 

First we apply Weyl's theorem with C = A and D = xB to conclude that 

{k: X k (A + xB) G I for some x G J} C {k : A*(A)e/}. (50) 

Suppose now that A^(A + jco#) G I for some value xq. Then using Weyl's theorem 
with C = A + xqB and Z) = (x — xq)B we obtain that 

X k (A + xB) ^7 for \x-x \ > a 1 (51) 

Combining (|50l) and (|5TI) we obtain the desired bound (|44l) . 

□ 



4 Renormalization of tadpole contributions 

The aim of this section is to set up the appropriate resolvent expansion that will 
be used in the proofs of Lemmas Q] and [21 namely to obtain decompositions (fT9l) 
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and (1241) . In particular, based on some combinatorial observation, equation (|60l ), 
renormalization of tadpoles is done in both real (case ^Y) and momentum (case 
ff) space. The estimates that control various terms in the resulting decompositions 
are established in Section [51 

We decompose as 

H^=H r + V, H r :=--A-a(p,E + ie), V := XVa + oip^E + ie) , 

where o(p,E + ie) is a solution of (TTOl) for {&) case. 
Respectively, for {JV ) case we decompose 

H^=H r + V, H r :=~A-Z, V:=XV a +Z, 

where E is a periodic extension of sigma defined in (fl4l) . Let 

R r ■= (H r - E - ie) r 1 . (52) 
We can expand 7? (defined in ©) into (truncated) resolvent series 

i? = Y,{- R rVyR r + (-R r V) N+l R. (53) 

7=0 

To handle the renormalization of tadpole contributions properly, we decide at 
which value of n to halt the expansion in (1531) individually for each contribution 
according to the following rule (to which we will refer as a stopping rule): If we 
open the brackets in (|53l , we obtain terms of the form 

R r QR r Q . . .R r 0R r 

where 6 is either — XV m , or — a(p, E + ie) I — E (whenever takes the later value 
we will refer to it as a bullet). Since o(p,E + ie) — 0(X 2 ), E = 0(X 2 ) for all 
permissible values of E, see Appendix|Bj one can unambiguously define the order 
I (in powers of A) of the particular contribution 

R r 0R r . . .R r QR§ , 

(with i?j being either R r or R) according to the following rule: Each factor of o 
counts as 2, while appearance of the random potential counts as 1, and we add 
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up all the exponents to get the order of the term. For instance, the order of the 
expression 

R r oR r XV m R r oR 

is 5. To illustrate this procedure we write down the expansion obtained in a case 
of N = 2: 

R = R r - R r aR - {XRrVtoR} = 

R r — R r oR — XR r VoyR r 

+ XR,y a R r oR + X 2 R r V C0 R r V (0 R, 

where the term in the curled brackets is the one we expanded according to the 
stopping rule. Note that the penultimate term is of order 3. It is not difficult to see 
(see Lemma 3.1 in |[6l for the proof) that for a general N we get 

Lemma 4. For any integer N we have a decomposition 

N-l N-l 

R = Y, A 'i R r+A' N R + B N R = £A, + A n R, (54) 

1=0 1=0 

where A ; = I, A\ is a summation over all possible terms of the type 

R r 8R r e...R r (55) 
which are of the order I > 0, while 

B N = -A N -iO. (56) 
The quantities Ai andAjy are defined as 

A, = A'jR r , A N = A' N + B N . 

In order to explain the renormalization, we borrow the following paragraph 
from Section 3.1. of [6] for notations. 

For an integer let T# be a set {l,...,N,N + 2,...,2N+ 1}. Let II = IL7 
be a set of partitions of into disjoint subsets Sj of cardinality |5/| G 2N. Two 
partitions % = {Sj}J =l , %' = {S'j\J =x are equivalent, % = it', if they coincide up 
to permutation. For S C Tyy, let 

5 ^)= i nv*h (57) 
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where 8 X , x G Z is Kronecker delta function, and xs denotes the collection of 
{xi , i G 5}. One has an identity (see e.g. J51 Section 3.1 for details) 



E 



n ^ 



N 



m=l 7C={Sj}J =i 7=1 



(58) 



where C2/ < (c/) 2/+1 and C2 = Eco 2 = 1, provided assumption holds. The set 
Sj in the partitions % G IT can be of the special type: If 



Sj = {i,i+l} 



(59) 



we will refer to it as a tadpole, or a gate set. 

Now let %? denote a collection of disjoint sets {Sj} such that any Sj G 7lf is a 
tadpole, and the cardinality of it? is fc. Then any partition n can be decomposed 
as K = %l U {5} for some < k < N, where S satisfies (U SjeJt cSj) US = fl^. Note 
that we didn't require 5 to be a tadpole free set. We will denote by Kq a partition 
of Tm such that no Sj G 71q is a tadpole. Lemmas |5] and [6] below hinge on the 
following observation: 



N 



k=0 



E 



7ren: 



/e5 



Tien-. 



Note that the summation on the right hand side runs over the tadpole-free parti- 
tions. To verify (1601) one just need to make a straightforward check that all tadpole 
contributions on the left hand side cancel out exactly. 

To see why the renormalization of tadpole's contribution in conjuncture with 
the above stopping procedure is useful in both (&) and (jV) cases, we consider 
two different tracks for each one of them: 



4.1 Case (JY) 

Let P x denote the projection onto the set — x for x G kT? , where 0, k are intro- 
duced in Assumption (jV). Then we can partition the identity operator as 

/ = £ Px, 
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and P x P y = for x ^ y. Instead of estimating the matrix elements of the corre- 
sponding terms in expansion (1541 ) directly, we will consider norms of operators 
P x RP y for x,y E kl? . Clearly, if x' E Range P x and y 1 E Range P y , then 

\R{x',y')\ < \\P x RP y \\. 

To evaluate P x RP y we use Lemma |H We insert the partitions of identity between 
each factor of the resolvent in (1551 ), with the net result 

£ P Xo R r P Xi e(x l )P Xl R r P X2 e(x 2 )...P Xn _ 1 RrPx n e{x n )P Xn RrP Xn ^ (61) 

6,Xj<EZ 3 ; 
j=l,...,n 

where 0(x) is either — X(0 X ^, or — E (defined in Eqs. QT|and[T2l respectively). 
The index n here depends on the particular contribution in A',, but the order of all 
contributions is /. 

To estimate the typical size of P x AiP y we consider the matrix 

< y := E {P x A[P y ■ PyA\P x } . (62) 

The key technical lemmas are the following assertions: 

Lemma 5. We have 

^en,: Xi ekZ 3 -. Sjen 
n=no J jeYi 

l-l 2/+1 

x X[{P Xi R r P Xi+l @} R r P y R* n {@P Xi R* r P Xi+1 } , (63) 
i=0 i=l+2 

where we are using convention xq = x; X21+2 = J- 
Proof. 

We first observe that by definition of @, E and (fl4l) we have 

X 2 ®P x R r P x 3) = P X L. 
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Using this identity in the definition of A/, we can represent 

N 



P x AlPyPyA]P x = ^£(-1)* £ £ 



i£5 



n <w 



x n{^A^ +1 ^}^x ri . (64) 



i=0 



i=l+2 



Computing the expected value of the left and right hand sides with respect to 
randomness and using (|60l ), we obtain (|63l ). 

□ 



4.2 Case {&) 

The counterpart of Lemma[5]in this case is the following generalization of Lemma 
3.2 of (6l (applicable for the non correlated randomness, i.e. u(p) = 1). In what 
follows, we will use the short hand notation E(p) in place of e(p) — E — ie — 
o(p,E + ie), and E* (p) for the hermitian conjugate of the multiplication operator 
E(p) . The renormalized propagator R r in this case will be given by its kernel 



i2ii(z-w)p d P 



E(p) 



(65) 



The following assertion holds: 

Lemma 6. For A/ defined in Lemma El the function E |A/ (x, y) | 2 is a function of 
the variable x — y. Let 



*flAx-y) := E|A z (x,v)| : 



(66) 



then we have 
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X 



Ja d Pl+l d P2l+2 pr dpj 



£(^+2) Efo) j±l\ 2 E*( Pj ) 

Y[u(pj-p j+ i) £ ri c l^| 5 [j^Pi-Pi+l h (67) 



K=TZq 



where 



a 



-i2n(pi+pi +2 )-(x-y). 
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Proof. 

Let Vj}y be a random potential of the form 

v£(x) = £ me~ m u(x-i). 

Then — > V a in the strong operator topology as 8 converges to 0. Similarly, 

we can define quantities , R 8 , R 8 , and A 8 by replacing V a with V 8 . One can 
readily check that 

R 5 (x,y)^R(x,y); R 5 r (x,y)^R r (x,y), A?(x,y) -> A,(x,y) 

in the limit 8 — > 0. The advantage of working with the regularized random poten- 
tial is due to the fact that it is summable and therefore admits Fourier transform. 
Namely, we have 

t£{p) = u( P )w s (p), 

where 

neZ 3 

Since by (fTOb we have 

A 2 / \u(p — q)\ 2 R r (q)dq = a(p,E + ie) , 

we can express \A 8 (x, y)\ 2 (analogously to (l64b ) as 

\A?(x,y)\ 2 = 

,21 f M dpi+i d P 2l+2 d Pj dpj _ 

W** E(p l+l )E%p 2l+2 ) }} lE ( Pj ) ^E^p,) ^ P J+ l) 

7C=7^U{S} 
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where /3 := 2k{ — (p\ + pi+2) (p/+i + P21+2) -y}- It follows from (|60l ) that 



E(-i)* E n s \Lpt-pt+i) E 

k=0 Ken,-. Si€7C% \ieSi 



PI ^«(Pi-Pi+l) 



e n%i 5 e^-P'+i ' 

:eIL Sjen \ieSj J 



Ken 

TZ=TtQ 

where stands for the convergence (with respect to 5) in the distributional 



sense. Therefore, taking the expected value on the both sides of (1681) as well as 
8 — > limit (where we use the smoothness of the integrand), we arrive to the 
expression that coincides with (1671 ), up to the prefactor e'P instead of e ia in the 
integrand. But the product of the delta functions allows to replace /3 with a (see 
Subsection |5.1| below), hence the result. 

□ 

5 Proof of Lemmas U and [2] 

Proof of Lemma [7] 

We observe that for any xf E Range P x and y 1 E Range P y 

ElA^x',/)! 2 < IK.IIi, (69) 

where srf x ,y is defined in (l62l) and || • ||i stands for the maximum absolute column 
sum norm. Indeed, 

\M^,y')\ 2 = {x'\P x Ai\y){y\A]Px\x') < (x'\P x AiPyA^Px\x') , 

hence 

E|A z (y,y)| 2 < E(x f \P x A l P y A}Px\x f ) < ik^ii < 
Therefore, using Lemma[5l we obtain that for such x' and / 

E\Mx', y ')\ 2 < x 2l \\n 21 E E n%i 5 K) 



^en,: Xj ekZ 3 : Sjen 

= 1/ 

/ 2/+1 



x niM^iiii n nw^iiii. ( 7 °) 

i=0 i=l+\ 
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with convention xi + \ =y,xo= X21+2 = x. 

Now we are in position to use Lemma[9]to bound the products of the norms on 
the right hand side of (1701) . for all |e| < k/2 defined in this lemma and all E < Eq 
with Eq be given by d8]) as 

Xl\\P Xi RrP Xi+l \\ x if ll^^illl ^ \®\ 2l+1 UGE*(x h X i+l ), (71) 
i=0 i=l+l i=0 

where Gs*(x,y) is defined in ( TT71) and E* = E*^ is given in the statement of 
LemmadJ We remind the reader that Ge* (x,y) is positive for any xje 1? . 
Plugging (TnT) into (1701) we get the estimate 



E|A/(x,y 



/ ../m2 



21+1 



71=710 j€T t 



;=0 



< A 2/ |0| 2/+2 ||^|| 2/ £ £ n c !5 7 -! 5 K) li^fe^i). (72) 

The latter expression, however, coincides (up to the factor |0| 2/+2 ||^|| 2Z ) with 
the corresponding term for the random potential of the form 

Vco(x) := £ (Qt, 

investigated in [|6]| (c.f. Eq. 3.15 there). As a result, the bound (|20l) follows from 
Lemma 1.1 of 0. 

To get (T2TI) note that it follows from Lemma |4] that 

A N = A N (H r -E -ie) - A N _{L. 

We therefore obtain 



E\A N (x,y)\ < £ {(E\A N (x,z)\ 2 ) l/2 -\(H r -E-ie)(z,y)\ 

+ (E\A N . 1 (x,z)\ 2 ) 1/2 • \L{z,y)\) 
< £ (E\A N (x,z)\ 2 f 2 + 2\®\X 2 \\D\\ 2 £ (E|A^!(x,z)| 2 ) 1/2 



zeZ 3 : zeZ 3 : 
\z-y\<l \z-y\e® 



25 



provided A is sufficiently small, and where in the last step we have used (1105I ). 

It now readily follows from (1201) that for A is sufficiently small, the right hand 
side of the above equation is bounded by 



, / a 2 X N/2 r- E 

C' |0| /+1 ||^||VW!^ [C ln 9 E* -j=J e" ~ 



□ 



with some generic constant C' . As a result, we have proved (12TI) . 
Further, ([221) follows from the fact that Ao(x,y) = R r (x,y) and 

|^r(jc,y)l < 

by Lemma|9l But the application of Lemma[7]fhen shows the validity of (l22l) . 
5.1 Feynman graphs 

At this point we have to introduce some additional notation: 

Definition 1. We consider products of delta functions with arguments that are 
linear combinations of the momenta {pi,P2, ■ ■ -,P2n+2}- Two products of such 
delta functions are called equivalent if they determine the same affine subspace of 
T 2n + 2 = {p h p 2 ,...,p 2n+2 }. 

One can obtain new delta functions from the given ones, by taking linear com- 
binations of their arguments. In particular, we can obtain identifications of mo- 
menta. 

Definition 2. The product 8(Y,j a jPj) of delta functions A n forces a new delta 
function, if £ 7 ajpj = is an identity in the affine subspace determined by A^. 

One can readily see that in the integrand of rhs of ([67]) one has a forced delta 
function 8{p\ — pi + \ + pi +2 — Pu+i), the fact used in Lemma[6] 

A[ e*(x — y) is conveniently interpreted in terms of the so called Feynman 
graphs (the pseudograph, to be precise, since loops and multiple edges are allowed 
here). The graph, associated with particular partition % of T W; „ is constructed 
according to the following rules (see Figure [TJ and [2]): We first draw two line 
segments, each containing n vertices (elements of T, hn ). The vertices are joined 
by directed edges (momentum lines) representing the corresponding momenta: 
.,p n+ i and p n+2 , . . . ,P2n+2- To each line pj we assign a propagator F(pj), 



26 



with some given function F, save momentum lines p\ and p n +2, which carry ad- 
ditional phases e -&npv{x-y) an( j e -i2x Pn+2 -{x-y) ^ reS pectively. For % = {Sj}J =l 
we identify all vertices in each subset Sj as the same vertex (in Figure [H the 
paired vertices are connected by dashed lines). Note that thanks to the existence 




Figure 1: Construction of the Feynman graph, part I, n = 4. The corresponding 
delta functions are 8(p\ -P2 + P3-P4), $(P4 -ps +P9 ~Pw), 8(P2~P3+P6- 
pj), and 8{pi — pg). The last delta corresponds to the tadpole. Note that the sum 
of all momenta in the above delta functions gives a forced delta function 8{p\ — 
P5+P6 —Pio)> hence we can introduce the dashed lines connecting vertices 1, 6,7, 
and 12, identifying them as a single vertex. 

of the forced delta function 8{p\ — p/ + i + pi + 2 — Pn+2)^ we can identify vertices 
1, /, / + 1, 21 as a single one, and therefore one can think about the closed graph 
(with special rules that apply for momentum lines p\ and p/ + 2> mentioned above). 
To summarize, the outcome of this construction is a directed closed graph, which 
is called the Feynman graph associated with the partition it. The momenta in 
the graph satisfy the Kirchhoff's first law, that is the total momenta entering into 
each internal vertex add up to zero (if arrow faces outward the vertex, we count 
its momentum with a minus sign). A tadpole corresponds to the so-called 0-loop, 
that is some (directed) line of the graph claims one vertex as its both endpoints. 
For a given Feynman graph G, one can choose a particularly useful expression 
for the product of delta functions A^. Choose any spanning tree of G which does 
not contain momentum lines p\,pi+2- The edges belonging to the spanning tree 
will be called the tree edges (momentum lines), and all the rest are the loop edges 
(since an addendum of any loop's momentum line creates a loop). Let us enu- 
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merate the tree variables as «i,...,M£, and loop variables as wi,...,w„, with say 
w\ = p\ , W2 = P1+2 (note that k + n = 21 + 2). The number k of the tree momenta 
coincides with the number of the delta functions in A n . 




Figure 2: Construction of the Feynman graph, part II: Identification of the vertices. 
The tadpole corresponds here to 0-loop. 

One can check (see e.g. flU) that the product of delta functions A K is equivalent 

to 

k i 

Y[8(ui- J^flyWy), (73) 

(=1 j=l 

with 

{±1 loop that contains u\ is created by adding w; to the spanning tree 
otherwise 

The choice of the sign depends on the mutual orientation of w, and Wj. 
Proof of Lemma \2\ 

We establish the exponential decay of £^i,e( x — y) m \ x ~ y\ from the following 
analytic argument. To this end we generalize the arguments used in Section 4. of 

to the momentum-dependent self energy o (defined in Eq.fTOl). 

We start with some additional notation: Let M denote a rectangle formed by 
the points 

{ -l/2±zVS; l/2±n/5}, 

where the parameter 8 = 8ff was introduced in (|23l) . For a unit vector e G T? , we 
will decompose T 3 3p = p- e©p i , where p L £ T 2 . In what follows we will use 
the norm || ■ || OT ^ defined as 

H/IU^ := max sup \f{p)\. (74) 

e p-eeM, p L eT 2 
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We note that for u satisfying © and 8 sufficiently small, one has 



- 4 Vd. (75) 

We will show that for a general value of /, 

**lfi{x) < Ml M -e-^ x \ J hE *(0), (76) 

where 

V) 2 '+ 2 e(Pl+i) +E* e( P 2i +2 ) +E* / e * e(/7 .) + £* 

x e n^i 5 [j^Pi-Pi+i > ( 7v ) 

K=71q 

with a defined in Lemma [6] 

The expression £^i,e*(0) has been studied in Section 4 of 10. It was shown 
there that 

4,e*(0) < W-E* (c\n 9 (E*) -±=) • (78) 



Combining (1761) . (1781) and (1751) . we obtain (1251) . Also, since Aq was defined to 
be equal R r in Lemma[2l we get ^,£*(0) = (G.e* (0,0)) 2 . Equation ([27]) follows 
from the bound (f8~7T) on the free Green function. 

The relation (l26l) is obtained analogously to the derivation of (|2TI) in the proof 
of Lemma [TJ 

To prove that (l76l) holds true let us choose for any given x El? the index 
ye {1,2,3} such that 

Ixvl = max |jcj| . (79) 
fe{l,2,3} 

Then \xy\ > \x\/y/3. In order to extract the exponential decay of ^^(x) we first 
perform the integration in the right hand side of (|67l ) over the tree momenta, using 

Let us use the shorthand notation f° r a sum over ai l possible partitions in 
(1671) . c n for a product of the corresponding c$j, r n for the number of the delta func- 
tions containing the loop momentum w\ in the 7T's partition, and s n will denote 
the number of u terms involving w\ after the integration of tree momentums. 
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Let E(p) = e(p) —E — ie — a(p,E + ie). We have 



/2n+2 i 2n 

c-^ xn d Wt [I ^TTT EI Mi), (80) 

where E*(p) stands for either E(p) or E*(p), 4>' is a set of indices of loop mo- 
mentum that does not include w\, and qu Qj are some linear combinations of the 
loop variables in 4>'. Note now that 

= J dwf e - i27t ( w V x -i w l-e r )x r ) 

Without loss of generality, let us assume that Xy > 0. The integrand as a func- 
tion of w\ ■ ey is 1-periodic, analytic inside the rectangle := C_ r\M for suffi- 
ciently small E*. Moreover, we have 

Ree(p-iVdey) > e(p)-2n 2 8 (82) 

uniformly in p G T 3 , provided £ is sufficiently small, where we have used the 
definition © and 

s'm(a + ib) = sinacosh^ + z'cosasinhZ?. (83) 
Combining this bound with (11011) , we get 

min (\E(p)\,\E(p-iV5e Y )\\ > e(p) + E* ,p G T 3 . (84) 

Moreover, the periodicity of the integrand implies that the integrals over the ver- 
tical segments of M- coincide: 



J -i/2 i = \E»(wi+qi) fJi 



L ^^Um^^'^U^^ <85) 
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Therefore 



/ d(wi ■ ey) ]1 W( - S e- i2 ™r("ve r ) Yl a( Wl + Qj) 



JT-iVd -Li E»(wi +qi) f = \ 

r T K 

I d(wi-e Y )H 

- 71 i=i 



E (wi+qj — ie r \/8 



< 



\qs« e -Wy/*fl f diw.-e^fl 1 (86) 

Y ' JT fjt e{w\ + qi) + E* 



where in the last step we have used (f8~4l) . Using the estimate (1 1 1 1> again, we also 
have 

\E(p)\ > e( Wl +qi)+E*, p E T 3 . 

Putting everything together on the right hand side of ((80l) . we get the bound 
(ESI). 

□ 



A Bounds on the free Green function 

The free Green function Gg(jc,)>) was defined in (TT7T ). We have 
Lemma 7. Define the function \j/ a e hC^+) as 

y a (r) = g-^max^-^)^- 2 )/ 2 , (l+r)^ 
TTzen /or J > 3 anJ — 1 < £ < we /iave 

o < G £ (*,y) < Q 

/or a// xj£ Z . 



(87) 



Remark A.l. A similar statement is known to hold on Mr, K25\l . We are not aware 
of its lattice version in the existing literature. The positivity ofG(x,y) on the lattice 
is well known, so it is an upper bound we are after here. 
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Proof. 

In what follows, we will use the following properties of the function i/a for d > 3 
and E* sufficiently small: 



(a) 
(b) 

(c) 
(d) 



IVolU = 1; 



£ llf a (\x-y\) = ^foranyyeZ^; 



Vo(k±*|) < C(0)v/ a (r)forO<£<diam(0); 



2n+l rTjj 2n+l 

II Va(\xi-i-Xi\) < e -l^ + i--ol^__ ]"[ ^(l^i-^l 
i=l za ;=l 

Suppressing the subscript E in the free Green function, we have 

J2n(x-y)p dd P = f J2x(x-y)p_d*P_ 



G{x,y) = / e lzn V- y >P . / = / e 



e(p)—E Jjd e(p)—E 
Let w = x — y. For any given w G Z d let us choose y e { 1 , . . . , J} so that 

|w-e r |= max I w - e,- 1 . (88) 

ie{l,...,d} 

Then 

|w-e y | > (89) 



Note that 



[ dp— 4 e- i2np - w = [ dqe- i2nqw d(p-e r )—4 e -V*fr'i»*r) 

J e(p)—E J J 1/2 e(p)—E 

(90) 

where q stands for the d — 1 dimensional vector obtained from p by removing its 
7 component (for d = 1 the argument below becomes completely straightforward, 
so we will only consider d > 2). Without loss of generality, let us assume that 
w ■ ey > 0. Let 

e(q) = 2 £ sm z (Kp-e a ). 
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It is easy to check that the integrand as a function of p- e r is 1 -periodic, analytic 
inside the rectangle formed by the points 



{-1/2; -1/2 W%^; + 1/2} 



6</ 



6J 



for a sufficiently small value of — E > 0: Indeed, using sin(a + z7?) = sinacosh/> + 
z cos ashing one can check that for any —1 < E < and e satisfying 



0<£< 



e(q) -E 



6d 



we have 

Ree(p + iee r ) -E > (e(p)-E)/2, 

uniformly in q. Moreover, the periodicity implies that the integrals over the verti- 
cal segments coincide: 



-1/2+// 



e(q)-E 
bd 



-1/2 

Therefore 

'1/2 



d(p ■ e 



1 



YJ e(p)-E 



— i2K{p-eyw-eyj 



bd 



1/2 



d(p-e Y )—4 e -iMp-*m) . (91) 

e{p)-E 



[ d(p-e Y )—^—e- i2 ^ w - e y 
J -1/2 e(p)-E 



bd 



— l/2+i 



6rf 



e{p)-E 



-iLn{p-eyW-ey 



< 2e e r\^r 



jq)-E fl/2 



-1/2 



d(p ■ e 



V e(p)-E 



y/e(q)-E 



(92) 



where in the last step we used (1891) . We can consequently estimate the right hand 
side of © by 

f uVfSEI 1 
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To estimate the latter integral, we split T d 1 into B := {q G T d 1 : e(q) < —E} 
and~B:=T d - l \B. Then 



B 



i i V e (i)- E 1 f i i V=£ 1 



y/e(q)-E Jb ^E 

< QrHT(-£)M/2 ] (93) 

and 



s y/e{q)-E J~B ^e{q) 



< dqe 
J~b 



2\/2q 



2 



<c,«W'f ('-2>M-g (94) 

for d >2, and where in the penultimate step we have used Jordan's inequality. 
Summing up d93]) and (|94i we arrive to (f87l) . 

□ 

Another useful property of the free Green function is captured by the following 
assertion: 

Lemma 8. For all E < and allT? 3 x^O we have 

1 < Ge(0,x) < 6 _ 2jBi (95) 



6-2£ G £ (0,jc + e) 
where e £ 7? is any unit vector. 
Proof. 

Note that for x ^ one has 

(0\(- 1 -A-E)(~A-E)- 1 \ X ) = 0. 

Inserting the partition of identity / = L z ez 3 \ z )( z \ between the operators on the 
right hand side, and using ([2]), we obtain the well known identity 

£ G E (-e,x) = (6-2E)G E (0,x). 

eeZ 3 : |e|=l 
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By translation invariance, Ge(— e,x) = G£(0,x + e). Since Ge(0,x) > for any 
x G Z 3 , we arrive to (1951) . 

□ 



B Properties of the self energy o 



B.l Properties of the solution of (TTOb 

In this section we establish the existence, periodicity, and analyticity of the self 
energy operator o(p, E) introduced in (ITOT ). We will use the following inequalities 
(that can be deduced from lfT4l ): 



Jt 3 e(q 



<2 



1 



"(//) -Iv 
To prove the existence, we introduce the space 



d 3 q- 

; (e(q) + £ 2 ) 2 £ 



1 

< - 



(96) 



L(T 3 ) = {/ : T 3 -> 



< 



oo , / is real analytic} . 



(97) 



We define a map T £ : L(T 3 ) — > L(T 3 ) pointwise as 

(T E f)(p)=^[ d\ { \ a{p - q) \\, r . 

We have T e Bp(0) C Br(0), where 5/3(0) is a ball (in || ■ ||oo topology) of radius /3 
centered at the origin, and 

J8 = 2A 2 ||«||i. 

Indeed, for/ €5^(0), 



Jt 3 



i«(p-^)r 



e(q)-E-i£-f(q) 



<X 2 f d\^<2X 2 \\u\\l. 
Jt 3 e(<?) 



Consider now the ball -By(O) of the radius 

7 := min(-£-2A 4 ||ii||i,2A 2 ||iC 
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Then T £ is a contraction on By(Q). Indeed, let f,g G By(0), then 



d q 



l"(p-?)| 2 |/(?) 



i= • \e(q)-E-ie-f(q)\\e(q)-E-ie-g(q)\ 

1 



< A 2 !!/-^^ / d 3 q 
Jt 3 



c 2 



< 



(e(q) + 2C*V) 2 " ^ 



= 11/- 



(98) 



where we have used (1961 ) in the last step. Hence by the Banach fixed point theo- 
rem, the self consistent equation (flOl) has a single valued solution o(p,E + ie) for 
all p E T 3 and all 

E < E := -2X 2 \\u\\l-2X 4 \\u\\i. 



The function o(p,E + ie) satisfies 

HctIL < min(-£-2A 4 ||i2||i,2A 2 ||w|| 2 ' 



(99) 



Next we establish 1 -periodicity of the above solution (in the real space). To 
this end, we note that since u(p — q) = L„ e z 3 u{n)e lm ^ p ~ q ' n , we have 

\u(p-q)\ 2 = £ u{m)u{n)e l7Zi{p - q){n - m) . 

Hence 



T 3 



d q 



\u{p~q)Y 



e(q)-E-i£-f(q) 

--X 2 V u(m)u(n)e 27tip{n - m) [ d 3 q 



2niq(m—n) 



e(q)-E-i£-f(q) 



(100) 



and periodicity of o(p,E + ie) follows from the periodicity of e 27tl P( n - m ) . 

Finally, we show analyticity. Fix a unit vector e G Z , and let p e := p ■ e, 
n e :=n-e for n GZ 3 . Using d96l) one can readily check that 



d k a(p,E + ie) 



4 

Pe 



< 2A 2 £ u{m)u{n){2n) k {n e — m e 

m,n£.1? 



< CX 2 A- k+3 k\. 



for all p e G T, and where A is given by ©, and C is some generic constant. This 
implies that o(p,E + ie) is real analytic in p e variable and admits the complex 
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analytic continuation to the rectangle M introduced in the paragraph followed by 
(|74l ). It follows from (flOl) and (|96l ) that we also have in this energy interval the 
bound 

||<7||oc,^ < 2A 2 ||w|| <xv ^, (101) 
with the norm || • \ x m defined in (|74|) . 



B.2 Properties of the solution of (1141) 



We proceed as in the previous subsection. We will be interested in the range of 
energies satisfying 

E < -k, (102) 

with K = An X 2 ||w|| 2 . It follows from the block diagonal structure of the operator 
E defined in (fT2~l) that for any pair Ei, E2 of such matrices (which correspond to 
<7i, 02, accordingly) we have ||Ei — E2II = ||<7i — 02 1 1 - Consider a ball 



B := {o eM nn : ||ct|| < k/2}. 



T : M n „ -+M n „, To := X 2 DSD, 



and a map 



where S is defined in (1131) and D in the paragraph followed by (fT2l) . We claim that 
TB C B. Indeed, for any (TGfi, we have 

||r<r|| < A 2 ||D|| 2 ||5||. 

To estimate ||5||, we observe that by construction of S, 

-1 



P(-A/2-E-i£-Zy l P 



where P is a projector onto supp 0. For the energies E that satisfy (11021 ) and o G B 
we have 

K(L) := Re(-A/2 -E-ie- E) > -A/2 + k/2. (103) 

But for an operator = A + z5 with positive A and hermitian 5, and a hermitian 
operator i 7 we have 



\F(A + iB)- l F\\ = FA- l l 2 (I + iA- l l 2 BA- l l 2 y l A- l l 2 F 



< 



FA- 1 / 2 



(I + iA-^BA' 1 / 2 )- 1 



A -m F 



< Wfa^fI 
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Hence 

P(-A/2-E-i£ 



< \\P (Re(-A/2-E-ie-Z))~ l P 



< 



p{-A/2+K/iy l p 



< 2n. 



where in the last step we have used the fact that the norm of the matrix is dom- 
inated by its trace norm, the positivity of P(— A/2 + k/2)~ 1 P, and bound ( |96| ). 
Putting everything together, we get 

||r<r|| < 2nX 2 \\D\\ 2 = k/2. 

By Brouwer's fixed point theorem the map T then has at least one fixed point in 
B. Since we are interested in proving the uniqueness, we will show that T is also 
a contraction on B. To this end, let 0\2 £ B. Then using the second resolvent 
identity 



- A/2 -E-ie-L 2 y l P 



(104) 



\Toi-To 2 \\ < A 2 ||D|| 2 ||Ei-E 2 || x 
P (-A/2-E-ie-Ei) -1 
Now observe that 
P(-A/2-E-ie-T.iY l (-A/2-E + ie-Ziy 1 P 

= PK^^iL^il + iB)- 1 K-\Li) (I-iB)- l K- l l 2 (U)P, 

where B is a self adjoint operator. Using (11031) . we can bound the right hand side 
(in the operator sense) by 

-PK- l l 2 {T.{){I + iB)- A (7-iB) _1 ^- 1/2 (Ei)P < -PK-\JL{)P < -, 

jA K K 

where in the last step we have used (|96l) . As a result, we have obtained the bound 



P (-A/2-E-ie-Li] 



< 



1 



K nX 2 \\u\\1' 
Using it and its analogue for E2 in (11041) . we arrive to 

1„ 

TCi — TO2W — _ O'l — °2 < °1 — °2 II) 
n 

hence T is a contraction on B. In summary, we have shown that ([14]) has a unique 
solution in the above energy interval, and 

||E|| = ||ct|| < 2|0|A 2 ||D|| 2 . (105) 
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B.3 Properties of R r in {jV) case 

Here we will consider the properties of the Green function R r (x,y) defined in (1521) 
where E satisfies ( TT5T ). The following assertion holds: 

Lemma 9. Let Eq be given by ®. Then for |e| < k/2, and all E < Eq we have 
\Rr(x,y)\ < (x\(-A/2-E+E Q /2)- l \y). (106) 

Proof. 

We first expand R r in Neumann series 

oo 

R r = G£((E + ie)Gy, 
;'=o 

with 

G := (-A/2-E)" 1 . 

Since ||E|| < k/2, the series converges absolutely for E < Eq and |e| < k/2. To 
estimate each individual term in the expansion, we insert partitions of identity 
I — £ ze ^3 \z) (z\ after each operator in the product. We obtain 

(x\G(ZGy\y) = £ G(x,zi)imz 2 i-i,Z2i)G(z2i,Z2i+i), (107) 

with a convention Z2j+\ = y- It follows from the construction of E that 

E(;c,;y) = Oforjc — y^®. 

Also, by (fT5l) we have 

|^(x } y)| < k/2. 

Using these bound together with the estimate (1951) to estimate the left hand side of 
(1107b . we get 

|<*iG(E*y»| 

< (|{(6-2 J B)«|0| + l}) J £ G(x,Z!) nG(z2/,Z2/+i) 

= (x\& +1 \y). (108) 
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Hence 



\G r (x,y)\ < (x\GY,(-E /2) j G j \y)= (x\ G (I + Eq/2 ■ G)~ ! \y) . 
i=o 

But 

G(I + E Q /2-GY l = (-A/2-E+E /2)-\ (109) 
hence the result follows. 

□ 



B.4 Dipole single site potential 

Here we consider a special case of the single site potential u c i defined in (fT6l) Let 
T £ be the same map as the one defined in (|97| ). Then \u(p) | 2 = 4sin 2 (7Tp • e\), and 
for the even function / we have 



JT 3 e(q)—E — i£ — 



e(q)-E-ie-f(q) 
+ 4A 2 / d\ , sin2 (^- g (110) 

We will consider the energies E that satisfy 

£ < -(1+A)A 2 . 

The subspace G of L°°(T 3 ) consisting of the functions f{p) = A + Bsin 2 (np ■ 
e\) is clearly invariant under the map T £ . Let G' denote an open subset of G' 
characterized by \A\ < A 2 , \B\ < 14A 2 . We then have T E G' c G'. Indeed, for 
/ G G ; , we can estimate the two terms on the right hand side of (II 101) using two 
bounds below, that hold for A sufficiently small: 

3 | cos(2tt<7 - ei) I f 3 1 +2sin 2 (^p-ei) 



V ^\e(q)-E-ie-f(q)\ Jv (1-|B|)«( 4 ) 

where in the penultimate step we have used the symmetry of the integral with 
respect to spatial directions {1,2,3} and in the last step we have used (|96l ). The 
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second estimate we need is 

3 2sin 2 (7z;g-ei) f 3 2 sin 2 (7T<5r ■ e\) 1 1 



T 3 d q \e( q )-E-ie-f(q)\ < (l-\B\)e(q) 3(1 - < 2" 

(112) 

Combining these two bounds we obtain 

|(r £ /)(p)| < 14A 2 sin 2 (^-ei)+A 2 , (113) 

hence T e / G G' . Since G' is a compact convex set, one can use Brouwer's fixed 
point theorem to conclude the existence of the fixed point (in fact one can use this 
technique to show the existence of the fixed point for all negative values of E). 
However, we also need a uniqueness of the fixed point, so we proceed to prove 
that T £ is a contraction on G' . To this end, let us introduce a norm on G: 

\\f\\ G = \A\ + X\B\, fotf = A + Bsm 2 {%p-e{). 
Let / = A + Bsin 2 (7tp -ei),g = C + Dsm 2 (7tp ■ e\) G G', then the straightforward 



computation similar to the one done in (1981 ) gives 

\\(Tsf)(p)-(T E g)(p)\\ G 

2 f , (\B - D\ + \A-C\) sin 2 (Ttq-eC 



< 4A 2 / d 3 q 
+ 4A 3 ! d 3 q 



((1- 


•]B\)e(q) + V 


)((1- 


\D\)e(q)+V) 




\A-C\ + \B- 


D\ sin : 


l {llq-e\) 


((1- 


\B\)e(q)+V) 


((1- 


\D\)e(q)+V) 



< 20A 2 |5-D| +5A 3/2 |A-C| < 20k\\f-g\\ G , (114) 

for A small enough. We have used (l9~6l) in the penultimate step. Hence by Banach 
fixed point theorem, the self consistent equation (flOl) has a single valued solution 
o(p, E + ie) for all p e T 3 and all 

E < E d := -(1+A)A 2 . 

Since o G G' we have 

o(p,E + ie) = A + Bsm 2 (np-ei); \A\ < A 2 , \B\ < 14 A 2 . (115) 

It follows from the functional form of o(p,E + ie) that for any unit vector e G 1? 
the function a is 1-periodic, analytic in p e := p ■ e (in fact it is a constant unless 
e = e\). Let E* d be a parameter that satisfies < E% < E d — E and let 



(E d -E-E*)/2. 
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Then using (|83l) and (II 151 ) we deduce that for an arbitrary p E T 3 we have 

Re(e(p + i8e)-E-ie-o(p + ide 1 E + ie)) > (1 - 5A 2 ) (e(p) + E* d ) . (116) 
Proof of Proposition [7] 

Let A: [— L,L] 3 HZ 3 , A + : [-L-1,L+1] x [-L,L] 2 nZ 3 andlet£2 A := x fceA+ M. 
By the standard arguments (c.f. the discussion in Section 6 of lfT5l ) it suffices to 
find a configuration co E Q. a, for which mine (H£) < — 2X 2 + 0(X 3 ) + o(l). Here 
o(l) is taken with respect to the L variable. We choose CO in such a way that for 
x E A, V^(x) = —2/1 for x • e\ =0 and (x) = otherwise. Clearly, the bottom 
of the spectrum of converges, in the limit L — >■ oo, to inf o(H) , where the latter 
operator acts on the whole Z 3 as 



H = --+V, 



A 
2 
with 

V(x) = —2X for x ■ e\ = and V'(jt) =0 otherwise. 

Readily, inf <j(H) < mina(^), where H is a one dimensional restriction of H to 
the e\ direction. However, H is a rank one perturbation of the free Laplacian. It 
follows from the rank one perturbation theory that E m := mina(//) is given by 
the solution of the equation 

2 A = Goo(E m ) , 

where G is a free one dimensional Green function. Using the Fourier transform, 
the above equation can be rewritten as 

I f dq 
2 A ./• 2sin 2 (7tq)-E m 



Finally, since 

dq f°° dq 



+ 0(1) = —— + 0(1) 



h 2 sin 2 (7tq) - E m 7— 2k 2 q 2 - E m V~2E„ 

which holds for E m < 0, we obtain 

E m = -2X 2 + 0(X 3 ). 

The rest of the argument coincides with the one of Theorem Q] for the case. 
One just need to replace the subscript 6 with d everywhere in the proof. 

□ 
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